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In terms of spinless fermions and spin waves, we describe magnetic properties of a spin-^ 
ferromagnetic-antiferromagnetic bond-alternating chain which behaves as a Haldane-gap antifer- 
romagnet. On one hand, we employ the Jordan-Wigner transformation and treat the fermionic 
Hamiltonian within the Hartree-Fock approximation. On the other hand, we employ the Holstein- 
Primakoff transformation and modify the conventional spin-wave theory so as to restore the sub- 
lattice symmetry. We calculate the excitation gap, the specific heat, the magnetic susceptibility, 
magnetization curves, and the nuclear spin-lattice relaxation rate with varying bond alternation. 
These schemes are further applied to a bond-alternating tetramerized chain which behaves as a 
ferrimagnet. The fermionic language is particularly stressed as a useful tool to investigate one- 
dimensional spin-gapped antiferromagnets, while the bosonic one works better for ferrimagnets. 

PACS numbers: 75.10.Jm, 75.40.Cx, 75.40.Gb 
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I. INTRODUCTION 

Haldanei*^ sparked renewed interest in one- 
dimensional Heisenberg antiferromagnets, predicting 
that their low-energy structures should qualitatively 
vary according as the constituent spins are integral or 
fractional. A magnetic excitation gap immediately above 
the ground state, which is referred to as the Haldane 
gap, was indeed observed in quasi-one-dimensional 
spin-1 Heisenberg antiferromagnets such as CsNiCla 3 
and Ni(C 2 H 8 N 2 ) 2 N0 2 (C104). 4 ' 5 A rigorous example of 
such a massive phase was also given theoreticallyi^ii 
Significant numerical effort s^iiSiiiiiS were dovoted to 
detecting the Haldane gap in the higher-spin systems. 
Competition between massive and massless phases 
in low-dimensional quantum magnets was extensively 
studied especially by the nonlinear-sigma-model quan- 
turn field theory 13 i 14 i 15 i 16 i 17 i 18 i 19 i 2Q i 21 i 22 i 23 and a wide 
variety of spin gaps— energy gaps in magnetic excita- 
tion spectra— were further predicted. There followed 
stimulative findings, including quantized plateaux in 
zero-temperature magnetization curves ^S^L gap 
formation in coupled spin chains 2 ^^ and the dra- 
matic crossover from one- to two-dimensional quantum 
antiferromagnets^ and an antiferromagnetic excitation 
gap with a ferromagnetic backgroundi 3 ^ 3 ^ 3 ^ 3 . 

Besides the sigma-model study, analytic approaches 
played a crucial role in revealing the nature of Haldane- 
gap antiferromagnets. The valence-bond-solid modeli 3 ^^ 
stimulated considerable interest in matrix-product 

of the Haldane phase. 
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The Lieb-Schultz-Mattis theorem 4 ^ was generalize 
clarify a mechanism for gap formation in a magnetic 
field. However, these arguments were essentially re- 
stricted to the ground-state behavior and can hardly 
be extended to finite-temperature properties. Numerical 
tools such as quantum Monte Carlo and density-matrix 
renormalization-group techniques are indeed useful for 
such a purpose, but an analytic strategy is still indis- 



pensable to low-temperature thermodynamics especially 
of spin-gapped antiferromagnets, where grand canonical 
sampling is hardly feasible numerically. Then we are led 
to describe massive spin chains in terms of conventional 
languages such as the Jordan-Wigner fermions and the 
Holstein-Primakoff spin waves. 

The Jordan-Wigner transformation is an efficient ap- 
proach to low-dimensional quantum magnetism. Spin-i 
arrays with uniform 4 ^ and alternating^ 3 ^ 4 ^ 5 " antiferro- 
magnetic exchange interactions between nearest neigh- 
bors were thus investigated and their energy structures, 
magnetization curves, and thermodynamic properties 
were indeed revealed well. Two-leg antiferromagnetic 
spin ladders were also discussed within this schem o 46 ! 47 
and the interchain-coupling effect on the lowest-lying 
excitation was elucidated. More refined fermionization 
was further proposed for coupled spin chains. Ordering 
spinless fermions along a snake-like path, Dai and Su 4 ^ 
succeeded in interpreting massive and massless excita- 
tions with varying number of the ladder legs. Their idea 
was generalized to investigate zero-temperature magne- 
tization curves^ and thermodynamic quantities^ - In 
such circumstances, we consider fermionizing an spin-i 
ferromagnetic-antiferromagnetic bond-alternating chain, 
which converges to the spin-1 antiferromagnetic Heisen- 
berg chain as the ferromagnetic coupling tends to infinity 
and therefore reproduces many of observations common 
to Haldane-gap antiferromagnets. 

Bosonic theory has significantly been developed for 
one-dimensional quantum magnets in recent years. While 
the Schwinger-boson mean-field theory is unable to dis- 
tinguish fractional-spin chains from integral-spin ones, it 
is still useful in predicting the asymptotic dependence of 
the Haldane gap on spin quantum number— and explain- 
ing quantum phase transitions of Haldane-gap antiferro- 
magnets in a field. 52 The Schwinger-boson representa- 
tion was further applied to ferrimagnetic spin chains^ 3 * 5 ^ 
and ladders— It was a major breakthrough leading 
to the subsequent development of the spin-wave the- 
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ory in low dimensiono 56 i 57 i 58 i 59 i 60 i 61 i 62 that TakahashiS 
gave a spin-wave description of the one-dimensional fer- 
romagnetic thermodynamics introducing an additional 
constraint on the number of spin waves. This mod- 
ified spin-wave scheme was further applied to spin- 
gapped antiferromagnets^S^S^ and qualitatively im- 
proved for one-dimensional ferrimagnets^S&Si The anti- 
ferromagnetic modified spin-wave theory is less quan- 
titative than the ferrimagnetic version^MS but it en- 
lighteningly interpreted novel observations such as the 
temperature dependence of the Haldane gap 64 i 65 and 
the field dependence of the nuclear spin-lattice relax- 
ation rateS Such spin-wave understanding is well sup- 
ported by other analytic descriptions^SiSLS As for 
finite-temperature calculation of spin-gapped antiferro- 
magnets, the Schwinger-boson mean-field theory is of 
no use, while the modified spin-wave theory maintains 
its validity to a certain extent^* Thus, we apply the 
modified spin- wave scheme to the spin-i ferromagnetic- 
antifcrromagnetic bond-alternating chain with particular 
emphasis on a comparison between fermionic and bosonic 
descriptions of spin-gapped antiferromagnets. 

Our theoretical attempt is much motivated by 
existent bond-alternating chain compounds such 
as IPACUCI3 (IPA = isopropylammonium = 
(CH 3 ) 2 CHNH 3 )S and (4-BzpipdH)CuCl 3 (4-BzpipdH = 
4-benzylpiperidinium = C^HisN)^ These materials 
behave as spin-1 Haldane-gap antiferromagnets at low 



temperatures^ 



75.76,77.78 



while such spin-1 features are 



broken up into paramagnetic spin ^'s with increasing 



temperature 



79,80.,' 



Besides the thermal crossover from 



quantum spin l's to classical spin ^'s, their enriched 
ground-state properties^2i2iSii&2& and novel edge 
states^ 7 , are of great interest to both theoreticians and 
experimentalists. 



II. FORMALISM 

We consider the ferromagnetic-antiferromagnetic 
bond-dimeric spin-i Heisenberg chain, whose Hamilto- 
nian is given by 



n 



N 

^(•/AF'S^ri-l ' S2n ~ ^F^ri ' S2n+l) 



-9PbH (S^n-l + 5|„) 



(2.1) 



The ground-state properties^2i22i2Si2i and low-lying 
excitations 92 of this model were well investigated by nu- 
merical tools and variational schemes. In particular, 
the string order parameter originally defined for spin-1 
Heisenberg chains^ was generalized to this system£M2i22i 
and the breakdown of a hidden Zi x Zi symmetry was 
extensively arguedi^S As the ferromagnetic coupling 
tends to infinity, the string order remains finite and the 
Haldane gap converges to that originating in decoupled 
singlet dimers. 



On the other hand, the thermodynamic properties have 
much less been calculated so fa»2L2i and there is no guid- 
ing theory for extensive experimental findings. Employ- 
ing two different languages, we calculate various thermal 
quantities and give rigorous information on their low- 
temperature behavior. 



A. Fermionic Approach 

In accordance with the bond dimerization, we intro- 
duce two kinds of spinless fermions through the Jordan- 
Wigner transformation 



^2n = blexp 



m—l 

n-1 



irr Y a L a m + b ™ br ' 



^2n-l — a n a n 7^1 — b\jj n — . 



(2.2) 



Decomposing the fermionic Hamiltonian at the Hartree- 
Fock level, we obtain a mean-field Hamiltonian as 

Huf = Eq + (Jaf — Jf) 

N r 1 1 

X H { db ~ o) " " + ( da ~ 2)^' 

N r 1 

+Jaf^ (~ -pxe^a^bn +H.c. 

n— 1 L 

~ Jf 5Z (2 ~Pp) b n a n+l +H.C. 



N 



-gfi B H Y ( a i a n + blb n 



(2.3) 



where d a = (a^a^w, d b = (b^b^BF, Paf = (bl,a n )n F , 
p F = (aJ t+1 6„) H F, and 



En 



Jaf(\p 



af| 2 — d a d b 



-Jf{ \pf\ 2 - d a d b 



N, (2.4) 



with (• ■ -)hf denoting the thermal average over the 
Hartree-Fock eigenstates. Defining the Fourier transfor- 
mation as 



b n = -±=j2 eikin+mb *> 



(2.5) 



and then a unitary transformation as 



bk 



u k v k e" 



v k e' 



Oik 

1% 



(2.6) 
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where 



Uk = 



Ik 



2 v Vv 2 + \ik\ 2 




V 



JA F (^-PA F )e [k/2 



Pf e 



-ik/2 



£ = - {Jaf - Jf) {d a + db - l) , 



(2.7) 



and twice the lattice constant is set equal to unity, we 
can diagonalize the Hamiltonian as 



hhf = e + J2 (44°* + £ k Ptfc) . 



where the dispersion relations are given by 

4 = £±vV + M 2 -<?mb#. 



(2.8) 



(2.9) 



In terms of the fermion distribution functions n fc = 

[e e fc /' bT + the internal energy, the total magneti- 

zation, and the magnetic susceptibility are expressed as 



E — Eq 
X 



E E «. 

k <y=± 



(2.10) 
(2.11) 



k a=± 
\2 



respectively. Another quantity of wide interest is the nu- 
clear spin-lattice relaxation rate 1/T\. Considering the 
electronic-nuclear energy-conservation requirement, the 
Raman process usually plays a leading role in the relax- 
ation, which is formulated as 



1_ _ 47r(ff^ B ^7N) 2 
T i ~ KE m e- Em/kBT 



E 



-E m /k B T 



''|E(^™'^2ri-l + BnS^n) 

71 

X5(E m > - E rn ~ 



(2.13) 



where A n and B n are the dipolar coupling constants be- 
tween the nuclear and electronic spins, wn = InH is the 
Larmor frequency of the nuclei with 7n being the gyro- 
magnetic ratio, and the summation ^2 m is taken over all 
the electronic eigenstates |m) with energy E m . Assuming 
the Fourier components of the coupling constants to have 
little momentum dependence as J2 n e -A n = Ak — A 



and J2 n e B n = Bk — B, we obtain the fermionic ex- 
pression of the Raman relaxation rate as 

Y = 7r(gM ^ 2 7N)2 E + B2 + ubohQ* - e k) ] 

1 k,k' 

X ^n%(l-nl,)5(e%,- E %-hws). (2.14) 



<7 = ± 



B. Bosonic Approach 

Next we consider a single-component bosonic represen- 
tation of each spin variable at the cost of the rotational 
symmetry. We start from the Holstein-Primakoff trans- 
formation 



9+ 

°4n-4+r 

s. 



2S- 



O'T-.rii 



s 



a 



An—i+T 

s. 



4n-2+ 



T — S + b T:n b T:n , 



(2.15) 



where r = 1,2; that is, we assume the chain to consist 
of four sublattices. Under the large-S" treatment, the 
Hamiltonian can be expanded as 

H = -2( J F + Jaf)S 2 N + E x + E +H X +H + O^ 1 ), 

(2-16) 

where Ei and TLi give the 0(S l ) quantum corrections 
to the ground-state energy and the dispersion relations, 
respectively. The naivest diagonalization of the Hamilto- 
nian fHty . whether up to O^ 1 ) or up to 0(5°), results 
in diverging sublattice magnetizations even at zero tem- 
perature. In order to suppress the quantum as well as 
thermal divergence of the number of bosons, we consider 
minimizing the free energy constraining the sublattice 
magnetizations to be zerosS^LiS 



JV 

EE( 

71=1 T=l,2 



bl n b T:n ) = 4SN. (2.17) 



Within the conventional spin-wave theory, spins on one 
sublattice point predominantly up, while those on the 
other predominantly down. The condition (|2.17|) restores 
the sublattice symmetry. In order to enforce the con- 
straint H2.17fl . we first introduce a Lagrange multiplier 
and diagonalize 

n = n + Jaf^E E ( + 4:A:n)- ( 2 ' 18 ) 

ra=l t=1,2 

We define the Fourier transformation as 
1 x 



br:n 



-ifc(n-5/8+r/4) 



a T -.k, 



1 



V^„ifc(n-l/8+T/4U 
k 



(2.19) 
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^2:k 




-^3:fc 








~^3:k 








V bL 


J 




\ -1>Ih 




4:t 


$1* 


J 



and then the Bogoliubov transformation as 



Pl;k 

V pL J 

(2.20) 

where four times the lattice constant is set equal to unity. 
We determine the coefficients ipf k so as to diagonalize H 
up to the order of 0(5) and perturbationally take TLq into 
calculation^ Then the Hamiltonian (|2.16|) is written as 

Ex = -2J AF (1 + 7 + v)SN + JafJ2J2 w * > 

k <r=± 

E = 2 J AF [2A(A - 7 r) - (1 + l)A 2 - 1 r 2 - A 2 ] , 



Ho 



k T=l,2 

EE 

k r=l,2 



Jaf 

"K/Tirrel ~l~ ^quart j 

where 7 = Jp/ Jaf , 

^ = 5V(l + 7 + ^ 2 



(2.21) 



l + 7 2 + 2axfe, 

^ = [^- 7 r- ( i + 7 )4](i + ^)i±^ 
- 7 (4 + r)l^ + (4 ^)^ + 7Si ^ 



r = 



1 

4JV 



EE 



k a xkK 
°i \ 1 + 7 + " 



1 



fe <7=± 

(1 + 7 + v)o X kCos 2 \ ~ 7(7 + CTX^A^sin 2 ! 



Xfe + 7(7 + 2axfc)sin"f 



,1 



EE 



• 2 t 

7sm I 



47V 
1 



k a=± 



a = —YY 

AN L^i L^i 

k a=± 



°XkK 

1 07 \ 1 + 7 + ^ _ 



AS 



(2.22) 



with \l = ul/S and X k = [(1 + 7 + ^) 2 - sin 2 (fc/2)] 1 / 2 . 
Wirrel and Hquart in 7~to contain off-diagonal one-body 
terms such as Oi T -kl3 T :k and residual two-body interac- 
tions, respectively, both of which are neglected in the 
perturbational treatment. 

At finite temperatures we replace a\. k a T: k and f3 T . k f3 r :k 
by their canonical averages (a^.^r-.k) = n^-.k an< ^ 
\.ui respectively, which are expressed as 

[ e J AF (^+5^)/fc B T _ Here the La _ 

grange multiplier v is determined through 

(2.23) 



n f:k = 



k <y=± 



Then the internal energy and the magnetic susceptibility 
are given by2& 



— (j — a 



k <r=± 



3k B T 



EE"£K + 1 )> ( 2 - 24 ) 



k ct=± 



respectively, where E g — — 2( J F + Jaf)S 2 N + E\ + E . 



III. CALCULATIONS 

First we calculate the ground-state energy E g and the 
lowest excitation gap £g ap and compare them with nu- 
merical findings in Fig. ^ The spinlcss fermions are 
much better than the modified spin waves at describ- 
ing both quantities. As Jf goes to zero, the fermionic 
findings are refined and end up with the exact values 
E g /N = — 3Jaf/4 and E gap = Jaf- The modified spin 
waves considerably underestimate the spin gap. They 
can not distinguish massive spin chains from massless 
critical ones^ to begin with, but they are still useful in 
qualitatively investigating dependences of the Haldane 
gap on temperature and spin quantum number S 

Secondly we calculate the thermodynamic properties. 
Figure El shows the temperature dependences of the zero- 
field specific heat and magnetic susceptibility. Due to 
the significant underestimate of the spin gap, the modi- 
fied spin- wave description is much less quantitative than 
the fermionic one at low temperatures. Furthermore, the 
modified spin waves completely fail to reproduce the an- 
tiferromagnetic Schottky-type peak of the specific heat. 
Because of the Lagrange multiplier v, which turns out 
a monotonically increasing function of temperature, the 
dispersion relations (|2.22|) lead to endlessly increasing en- 
ergy and thus nonvanishing specific heat at high temper- 
atures. The spinless fermions succeed in reproducing the 
overall thermal behavior. The present approaches have 
the advantage of giving the low-temperature behavior an- 
alytically. Equation l|2.9|l shows that the dispersion rela- 
tion of the low-lying excitations reads 



4 - ±(£gap 

provided g^ B H < E gap , 



+ J A Fvk 2 ) 
where 



(3.1) 



2E, 



j afPaf + j fPf + 2J F JafPfPaf, 



gap 1 



JfPfPaf, 



(3.2) 



with pf — Repp — 1/2 and paf = RepAF — 1/2. Then 
the low-temperature properties are calculated as 



c r 
(gm) 2 N 



ttvJaf 



\knTJ 



knT 4 



Jaf _ e -E s 



irvk B T 



,/k B T 



(3.3) 
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These features are found in the antiferromagnetic Heisen- 
berg two-leg ladder as wel l 50 i 72 and can be regarded as 
common to spin-gapped antiferromagnets. The power- 
law prefactor to the activation-type temperature depen- 
dence, which can hardly be extracted from numerical 
findings, is essential in estimating the spin gap experi- 
mentally. 

Next we consider the total magnetization as a func- 
tion of an applied field and temperature. We compare 
the fermionic description of magnetization curves with 
numerical findings in Fig. [3] The spinlcss fcrmions again 
work very well. Quantum Monte Carlo sampling becomes 
less and less feasible with decreasing temperature, while 
we have no difficulty in calculating Eq. (|2.11(l even at 
zero temperature. The ground-state magnetization turns 
out to behave as M oc (H — He) 1 / 2 near the critical field 
g/iB-ffc = E gap £2" Magnetization plateaux of multi-leg 
spin ladders^, and mixed spin chains^ are also well inter- 
preted in terms of the spinless fermions. On the contrary, 
in the modified spin-wave theory, the number of sublat- 
tice bosons are kept constant and therefore we have no 
quantitative information on the uniform magnetization 
as well as the staggered one. Though the Schwinger- 
boson mean-field theory^i^SSiiffi which consists of a 
rotationally invariant bosonic representation, still works 
with an applied field and/or existent anisotropySiiSiiiS 
to a certain extent but rapidly loses its validity with in- 
creasing temperature^ 

Thus and thus, we are fully convinced that the spinlcss 
fcrmions are superior to the modified spin waves in inves- 
tigating quantum and thermal properties of spin-gapped 
antiferromagnets. Lastly in this section, we calculate the 
nuclear spin-lattice relaxation rate 1 jT\ in terms of the 
spinless fermions in an attempt to stimulate further ex- 
perimental interest in this system. If we again employ 
the approximate dispersion (|3.1|l at moderate fields and 
temperatures, k^T <C E gap — g^H, Eq. (|2.14|1 can be 
further calculated analytically as 



1_ ^ (g/i B ft7N) 2 (A + B) 2 e -E s ^/k B T 

9^bH ( Tiuj 
xcosh— — K 



V2fc R T/' 



(3.4) 



where Kq is the modified Bessel function of the second 
kind and behaves as Kq(x) ~ ln2 — 7 — \nx for < x <C 1 
with 7 being Eulcr's constant. Considering the signifi- 
cant difference between the electronic and nuclear energy 
scales (Tiwn ^ 10 _5 J), there usually holds the condition 
Huj-n <C k^T. At low temperatures, I/Xi also exhibits 
an increase of the activation type but with logarithmic 
correction, which is much weaker than the power correc- 
tion in the case of the susceptibility Such a pure spin- 
gap-activated temperature dependence of 1/Ti, which is 
shown in Fig. 0] should indeed be observed experimen- 
tally, unless magnetic impurities mask the intrinsic prop- 
erties. Equation l|3.4|l further reveals a unique field de- 
pendence of 1/Ti: With increasing field, 1/Tj first de- 
creases logarithmically and then increases exponentially, 



which is visualized in Fig. The initial logarithmic 
behavior comes from the Van Hove singularity peculiar 
to one-dimensional energy spectra and may arise from 
a nonlinear dispersion relation at the band bottom in 
more general. Therefore, besides spin-gapped antifer- 
romagnets, one-dimensional ferromagnets and ferrimag- 
nets may exhibit a similar field dependence&SiiSiiSi 
Relaxation-time measurements on spin-gapped chain an- 
tiferromagnets such as IPACuCLj and (4-BzpipdH)CuCl3 
are strongly encouraged. 



IV. BOND-ALTERNATING FERRIMAGNETIC 
CHAIN 



Before closing our comparative study, we briefly men- 
tion a bond-alternating but ferrimagnetic chain calcu- 
lated within the same schemes. We take another interest 
in the ferromagnetic-ferromagnetic-antiferromagnetic- 
antiferromagnetic bond-tetrameric spin-^ Heisenberg 
chain, whose Hamiltonian is given by 

AT 

4n-3 ' Sin-2 + S 4 n -2 ' $4n-l) 

71=1 

— Jp(Si n -l ■ Sin + S4n 1 SUn+l)] ■ (4-1) 
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Cu(3-Clpy)2(N 3 ) 2 (3-Clpy = 3-chloropyridine 
C 5 ClH 4 N)i2^ is well described by this Hamiltonia: 
and behaves as if it is a ferrimagnet of alternating spins | 
and I In the conventional spin- wave scheme, the spin 
deviations in each sublattice, (flJ.. n a T: ») and (6^. n & r: „), 
diverge in the antiferromagnetic ground state but stay 
finite in the ferrimagnetic one. Without quantum diver- 
gence of the sublattice magnetization, it is not necessary 
to diagonalize the effective Hamiltonian 12.18fl . In an 
attempt to keep the dispersion relations free from tem- 
perature, we may simply diagonalise the original Hamil- 
tonian (|4.1() and then introduce a Lagrange multiplier so 
as to minimize the free energy^ For ferrimagnets such 
an idea is much superior to the original antiferromagnetic 
modified spin- wave scheme^Si^LS 

Figure |S] shows the thus-modified spin-wave calcula- 
tions as well as the Hartree-Fock calculations in terms of 
the spinless fermions in comparison with numerical find- 
ings. The ferrimagnetic modified spin waves work very 
well, contrasting with the antiferromagnetic ones. They 
reproduce the Schottky-type peak of the specific heat 
and the ferrimagnetic minimum of the susceptibility- 
temperature product at intermediate temperatures. Al- 
though the modified spin-wave description of the anti- 
ferromagnetic increase of \F is somewhat moderate, it 
converges into the paramagnetic value S(S + l)/3 at 
high temperatures. The description is more and more 
refined with decreasing temperature and is expected to 
be accurate at sufficiently low temperatures & The T 1 / 2 - 
initial behavior of C and the T _2 -diverging behavior of 
X are both correctly reproduced^ Besides static proper- 
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ties, T\ measurementsiSi on a ferrimagnetic chain com- 
pound NiCu(C 7 H 6 N 2 6 )(H 2 0) 3 -2H 2 was elaborately 
interpreted in terms of the modified spin waves^£& 

On the other hand, the spinless fermions misread the 
low-temperature properties of ferrimagnetic chains. A 
fatally weak point of their description is the onset of 
a N eel-ordered state. With increasing Jp, the transi- 
tion temperature T c goes up and the applicability of the 
Hartree-Fock fermions is reduced. Indeed the fermionic 
description is not so bad away upward from T c , but it is 
much less complementary to numerical tools in ferrimag- 
netic systems. 

V. SUMMARY 

We have comparatively discussed fermionic and 
bosonic descriptions of the bond-dimcric Heisenberg 
chain as an example of spin-gapped antiferromagnets. 
The fermionic language is based on the Jordan- Wigner 
spinless fermions within the Hartree-Fock approxima- 
tion, while the bosonic formulation consists of con- 
straining the Holstein-Primakoff bosons to restore the 
sublattice symmetry. The spinless fermions well de- 
scribe both ground-state and finite-temperature prop- 
erties. The zero-field specific heat and magnetic sus- 
ceptibility behave as C oc (fc B T) -3 / 2 e~ f; £ £, p/ A:BT and 
X <x (fcBT)~ 1 / 2 e~ Bgap / fcBT , respectively, at sufficiently 
low temperatures, while the relaxation rate as \jT\ oc 
e-^p/ fcBT cosh(^B J ff/fc B T 1 ) [0.80908- \n{hcj N / k B T)} at 
moderately low temperatures and fields. On the other 
hand, the modified spin waves give much poorer findings. 



In particular, they significantly underestimate the spin 
gap and fail to reproduce the Schottky-type peak of the 
specific heat. The same schemes have further been ap- 
plied to the bond-tetrameric ferrimagnetic chain, where 
the modified spin waves work very well and are superior 
to the spinless fermions both qualitatively and quanti- 
tatively. The fermionic language is useful in describing 
disordered ground states and their excitations, whereas 
the bosonic one in depicting ordered ground states and 
their fluctuations. 

The modified spin-wave theory is fully applicable to 
higher-spin systems. The Jordan- Wigner transforma- 
tion can also be generalized to higher-spin systems?-^ 
where spin-1 chains, for instance, are mapped onto 
an extended t-J model of strongly correlated electrons. 
However, the double- graded Hubbard operators such as 
c n) f = (1 — c' n i Cn.jjcn,! demand that we should treat 
the fermion and boson degrees of freedom in the same 
footingiSiLiifliiiiiiiS The present naive fermionic repre- 
sentation is highly successful for spin-i gapped antiferro- 
magnets, including various bond-alternating and/or cou- 
pled chains. It is complementary to numerical tools espe- 
cially at low temperatures and allows us to readily infer 
both static and dynamic properties of spin-gapped anti- 
ferromagnets. 
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FIG. 1: The spinless-fermion (SF), modified-spin-wave (MSW), quantum Monte Carlo (QMC), and numcrical-diagonalization 
(Exact) calculations of the ground-state energy (the left) and the excitation gap immediately above the ground state (the right) 
for the bond-alternating dimerized chain, where L = 2N is the number of spins. 
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FIG. 2: The spinless-fermion (SF), modified-spin-wave (MSW), and quantum Monte Carlo (QMC) calculations of the specific 
heat (the upper three) and the magnetic susceptibility (the lower three) as functions of temperature for the bond-alternating 
dimerized chain, where L = 2N is the number of spins. 




g^ B ///7 AF g^H/J k¥ 8^bH/Jaf 



FIG. 3: The spinless-fermion (SF) and quantum Monte Carlo (QMC) calculations of magnetization curves for the bond- 
alternating dimerized chain, where L = 2N is the number of spins. 
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FIG. 4: The spinless-fermion (SF) calculations of the nuclear spin-lattice relaxation rate as a function of temperature (the left) 
and an applied magnetic field (the right) for the bond-alternating dimerized chain. 




FIG. 5: The spinless-fermion (SF), modified-spin-wave (MSW), and quantum Monte Carlo (QMC) calculations of the specific 
heat (the left) and the magnetic susceptibility (the right) as functions of temperature for the bond-alternating tetramerized 
chain, where L = AN is the number of spins. The Hartree-Fock fermions encounter a paramagnetic-to-Neel-ordered phase 
transition with decreasing temperature and the transition temperature is indicated by arrows. 



